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Abstract. - Motivated by recent experiments showing the compressive buckling of microtubules 
in cells, we study theoretically the mechanical response of, and force propagation along elastic 
filaments embedded in a non-linear elastic medium. We find that embedded microtubules buckle 
when their compressive load exceeds a critical value fc, and that the resulting deformation is 
restricted to a penetration depth that depends on both the non-linear material properties of the 
surrounding cytoskeleton, as well as the direct coupling of the microtubule to the cytoskeleton. 
The deformation amplitude depends on the applied load / > /c as (/ — fc)^^^. This work shows 
how the range of compressive force transmission by microtubules can be tens of microns and is 
governed by the mechanical coupling to the surrounding cytoskeleton. 
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The mechanical response of most eukaryotic cells de- 
pends on their cytoskeleton, a composite network of fila- 
mentous proteins [1]. Microtubules (MTs) are the stiffest 
of these cytoskeletal filaments, and they play an impor- 
t-ant role in organization of, and transport within the cell. 
Their mechanical rigidity allows them to support signif- 
icant stresses in the cytoplasm. These stresses can be 
highly inhomogeneous, with compressive/tensile forces di- 
rected along stiff MTs, permitting directed force transmis- 
sion and mechanical signaling over several microns within 
the cell. As with macroscopic elastic rods, however, even 
the comparatively rigid MTs cannot, on their own, with- 
stand as large compressive loads as tensile loads. This 
is because of the classical Euler buckling instability lim- 
iting the compressive force to a maximum value, which 
actually vanishes for long rods. It was recently shown, 
however, that even long MTs can bear large compressive 
loads, as a result of their coupling to the surrounding elas- 
tic matrix of the cytoskeleton [2] . This composite aspect of 
the cytoskeleton has important consequences for cell me- 
chanics and mechanotransduction [3-6] — the generation, 
transmission, and sensing of forces by the cell. 

Here, we develop a model for compressively loaded elas- 
tic filaments such as MTs embedded in an elastic con- 
tinuum. When their compressive load / exceeds a critical 
force fc, an oscillatory buckling of the filament is expected. 



with a wavelength depending on both the stiffness of the 
elastic filament and the shear modulus of the surrounding 
medium [2,8,9]. In the classical Euler buckling problem, 
and even in the presence of a (linear) elastic background, 
an elastic rod becomes unstable for f > fc- We include the 
non-linear elastic properties expected for the cytoskeleton, 
and show that the system is stable to supercritical loads, 
with a buckling amplitude that increases above thresh- 
old as 1/ — /c|^^^- In addition, both the buckling ampli- 
tude and the compressive load decay away from the point 
of force application in a way that depends sensitively on 
the longitudinal mechanical coupling of the filament to 
its surroundings. This suggests that the range of force 
transmission in the cell can be effectively controlled by 
microtubule-associating proteins that couple MTs to the 
rest of the cytoskeleton. Surprisingly, the experiments in 
Ref. [2] (see Fig. 6 therein) also found evidence that, de- 
spite the existence of a threshold force for buckling, the 
deformation of the microtubule was attenuated, suggest- 
ing a spatially varying force. 

In studying the mechanical response of intracellular 
MTs, we extend the classical buckling theory of rods in an 
elastic medium [2,8,9] to take into account the mechanical 
coupling of both longitudinal and transverse deformations 
of the MT to the surrounding nonlinear elastic cytoskele- 
tal network. We use a variational theory, numerical (con- 
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jugate gradient) minimization, and asymptotic analysis of 
the deformation of the MT. On applying a compressive 
load at one end, the microtubule bends and deforms the 
surrounding elastic network which resists this deforma- 
tion. These competing effects determine the characteris- 
tic force and length scales associated with the buckling of 
the MT. The nonlinear clastic properties and longitudi- 
nal coupling of the MT to its surroundings determine the 
attenuation of both the force and deformation of the MT 
along its arclcngth away from the point of loading. While 
the details of the deformation profile depend on the pre- 
cise form of the nonlinearity, the force penetration depth 
gcnerically decreases with increasing strength of the lon- 
gitudinal coupling between the MT and the medium. 

We consider a linear elastic theory for the MT and the 
surrounding medium and include the leading nonlinear- 
ity for the elasticity of the surrounding matrix. The MT 
is subjected to a compressive load / at one end denoted 
by a; = 0. The resulting elastic deformation energy can 
be written in terms of its transverse and longitudinal dis- 
placements, u{x) and v(x), respectively: 



-MO) 



-u (x) + —u[x) 
+ -u[x) 



dx 
dx. 



where v{x) = ^°°\u (y) dy as is required for an incom- 
pressible rod. The first term represents the axial com- 
pression energy released by bending. The second term 
corresponds to the bending energy of the rod, k is the 
bending rigidity. Classical Euler buckling, in which a 
long-wavelength bend of a rod of length L occurs at a 
threshold force of tt^k/L^, results from just the first two 
terms [8, 10]. The elastic constants a±, a|| correspond 
to the linear elastic transverse and longitudinal mechani- 
cal couplings between the MT and the surrounding actin 
network. The transverse coupling is directly related to 
the (linear) shear modulus G of the surrounding medium: 
a± = 47rG/ ln(A/a), where A is the characteristic wave- 
length of lateral displacement and a is a small wavelength 
cut-off. Adding just this third term leads to buckling on 
a short wavelength A = 2tt(!q, where io = (K/a^)^/* [2,8], 
and at a force threshold fc = Stt^k/X"^ = 2{na±y^^ that 
is significantly larger than for Euler buckling. The longi- 
tudinal coupling, in contrast, is not directly related to the 
elasticity of the surrounding medium, but also depends 
on the strength/degree of attachment between the MT 
and the surrounding network. Thus, it falls in the range 
< q:|| < 27rG/ln(A/a). Here, the lower limit corresponds 
to an uncoupled MT that slides freely along its axis, while 
the upper limit corresponds to the strong coupling of the 
MT to the medium along its length. Finally, the term 
proportional to /3 > represents the leading non-linear 
transverse coupling to the medium. In fact, both terms 
in the second integrand of Eq. (1) are non- linear in the 
transverse displacement u since the first of these is effec- 




Fig. 1: The deformation of a stiff elastic filament embedded in a 
nonlinearly elastic medium decays as an oscillatory exponential 
(from numerics). The material parameters are k, — P = a± = 
1.0, and ail = 0.1. 



tively fourth-order in u, due to the quadratic relationship 
between u and v. 

We examine the effect of the nonlinear terms in Eq. (1) 
by numerically determining the deformation of a MT in 
the surrounding matrix that is subjected to symmetric 
compressive loads applied at both ends: x = and x = L. 
In the limit of large L, the solution near x = is equiva- 
lent to that of a MT loaded at just one end, and the sym- 
metrized loading is simply for computational convenience. 
Using a conjugate gradient technique [11], we minimize 
the energy of the loaded system given by Eq. (1) with an 
additional term +fv{L) representing the load at .t = L 
and subject to hinged boundary conditions. Wc find that 
the buckling threshold fc and wavelength A = 2tt£o re- 
main approximately the same as in the linear system, but 
that the amplitude of the undulatory deformation u{x) 
just above threshold (i.e., for small / — /c > 0) is highly 
attenuated (see Fig. 1). The buckling amplitude scales 
with the applied force / as (/ — fcY^^ reminiscent of a 
supercritical pitchfork bifurcation [12]; it also scales with 
the nonlinear elastic constant /3 as a power law with ex- 
ponent — 1/2 as shown in the inset of Fig. 2. At a given 
force, the oscillatory deformation field decays with an ap- 
parently exponential envelope as one moves along the axis 
of the MT away from the loading point, with a penetra- 
tion depth or decay length that scales as {P/a^^Y^^ (Fig. 
3). 

Wc propose a variational ansatz for the form of the 
transverse deformation field based on the numerical so- 
lutions shown in Fig. 1: 



u{x) = uo exp (— .t/^) sin (gx). 



(2) 



Substituting this into the our energy expression, Eq. (1), 
we determine uq, i and q so as to minimize that en- 
ergy. For applied loads less than a threshold given ap- 
proximately by fc = 2(Kaj^)^/2, the minimum requires 
a straight MT: uq = q ^ 0. Above this threshold, the 
minimum is at a finite value of uq and q ~ l/^o- Once 
again, the buckling amplitude grows with applied load as 
uo ^ if ^ fcY^^ and decreases with an increasing nonlin- 
ear elastic constant f3 as as found numerically (see 
Fig. 2). At a given load above threshold, / > fc, the pen- 
etration depth of the oscillatory deformation field scales 
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Fig. 2: The amplitude uo as function of compressive load / 
for K = a± — 1.0, a|| = 0.1, ai P — 1 (circles) and /3 — 10.0 
(squares). Inset shows the amplitude as a function of /3, with 
f / fc = 1.25. In both the main figure and the inset, filled 
symbols correspond to the numerics, open symbols to the vari- 
ational calculation. The dashed line represents the prediction 
of the asymptotic analysis. 



as (/3/q!||)^/^ (Fig. 3). The results of the variational calcu- 
lation, therefore, corroborate all of the observations from 
the numerics. 

We nondimensionlize the energy in Eq. (1) by scaling x 
and all displacements w, u by .^0 — e.g., u{x) = u{x/£o)/£o. 
We also scale the energy by fo^o, where /o = {Ka±) 



U/2 



The resulting dimensionless energy £ 
ten as 



£ 

Mo 



can be writ- 



-#(0) 



1 // 



1 



2"" (*) + 2"(-*) 
-v{x) + -u(x) 



dx 
dx. 



(3) 



where 4> = f I fa is the dimensionless applied force and 
X = xjiQ. We have also introduced dimensionless ma- 
terial parameters e = q;||/q;j^ and 77 = 151^' ja^. Based 
on the estimates above for q;||,_l, we expect that e < 1/2. 
The local strain is approximately given by w' ~ u. For 
cytoskeletal networks, significant stiffening is expected for 
strains less than 1 [14,15]. Thus, we expect 77 ^ 1. Our 
calculations are all for 1)1 1 > 10, 0.5 < rj < 100, and 
0.001 < e < 0.1. 

In terms of the dimensionless variables, the ansatz in 
Eq. (2) becomes u{x) = exp (— i/A) sin (gi), with A = 
£/£o and q = loq. Using this to minimize Eq. (3), we 
find the dependence of A and ip on the applied load near 
threshold, (0 — (f>c)/(pc ^ 1, in the case of weak longitudi- 
nal coupling (e 4C 1), as follows. In this limit, A ^ 1 so 
we expand the energy Eq. (3) in 1 /A and keep the leading 
terms, which are of order A and eA'^. Dropping correc- 
tions that are suppressed by a factor of either e or 1/A^, 
we determine the minimum of the truncated, asymptotic 
expansion of the energy with respect to -0, A, and q. We 




Fig. 3; Decay length / as a function of the longitudinal coupling 
constant ay for k = P = a± — 1.0, ///c = 1.25. Inset shows 
the variation of the decay length with the nonlinear elastic 
constant /3 for k = a± = 1.0 and ay — 0.01 and f/fc = 1.25. 
Filled symbols in both the main figure and inset correspond to 
the numerics, open symbols to the variational calculation and 
the dashed lines to the asymptotic analysis. 



find q = q£Q = 1, and 



A 



1277 
e 



V' = ± 



3r] 



=)• 



(4) 



To lowest order in 1/A and e, the critical force fc is 
unchanged from that of the linear theory. Specifically, 
(j)c = 2 + 4A~^ In Fig. 4 we show results of the 

numerics (filled symbols), the variational analysis (open 
symbols) and the asymptotic analysis (lines) demonstrat- 
ing the excellent agreement of these approaches. 

The term in Eq. (3) represents only the initial non- 
linear correction to the elasticity of the matrix. We expect 
a stronger nonlinear response for crosslinked cytoskeletal 
networks. Prior theory and experiment suggest an inverse 
quadratic divergence of the stress versus strain [13-15] 
for F-actin networks. Motivated by this, we replace the 
transverse elastic terms in Eq. (4) a term of the form 
{u^jA)[\{u/us) + + \{u/us) -2]. This is chosen 

to reproduce the term in Eq. (3) to leading order. For 
small ft, in fact, this should be equivalent to Eq. (3) with 
rj = 2/{tj. With this more physical nonlinearity, we find 
a clear departure from the exponential envelope for buck- 
ling oscillations observed earlier. In Fig. 5 we show the 
amplitude of the buckled MT for Ug = 1.4; there appears 
to be a transition from a regime of slow amplitude decay 
near the point of loading to an exponential decay with 
the same rate as for the model in Eq. (3) with t] = l.The 
inset compares the actual decay envelopes for the sim- 
ple and higher order nonlinearities. The latter decreases 
the amplitude of buckling near the point of loading, how- 
ever, the envelope asymptotes to exponential decay with 
a decay length equal to that for the simple nonlinearity 
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Fig. 4: Collapsed data from Figs 2 and 3. The main figure 
shows the scaled amplitude as a function of the dimension- 
less force (f). The inset shows the dimensionless decay length 
A as a function of e/rj the relative strength of the nonlinear 
terms for longitudinal and transverse coupling. Filled symbols 
correspond to the numerics, open symbols to the variational 
calculation, and the dashed lines to the asymptotic analysis. 



of Eq. (3). The asymptotic approach to the exponential 
decay envelope (as predicted by Eq.(3)) at smaller defor- 
mation amplitude is expected given that the elastic energy 
shown in Eq.(3) reproduces the lowest order terms of the 
more physical elastic nonlinearities explored here. 

Finally, we note that we have assumed the MT to be 
incompressible. Finite compressibility introduces another 
decay length; the compressive strain on an elastic rod with 
compression modulus ^j. and longitudinal coupling to the 
network of a|| decays over a length scale yj^i[a\\ in the un- 
buckled limit. Treating a MT as an elastic rod of radius 
a ~ lOnm and Young's modulus IGPa [16], we estimate 
this decay length to be ^ 10//m for the largest a\\ ~ IkPa. 
At smaller e, where our our model is applicable to intra- 
cellular MTs. the actual decay length will be larger than 
this, likely spanning the whole cell. 

Motivated by observations of MT buckling in the cy- 
toskeleton, we have explored force propagation and com- 
pressive buckling of stiff rods embedded in model nonlinear 
elastic matrices. Consistent with prior experiments [2], we 
find oscillatory buckling with a decaying amplitude and fi- 
nite spatial extent I. We have focussed on the limit of weak 
longitudinal coupling q;||, where the penetration depth I 
is large compared with the buckling wavelength A. We 
estimate these lengths using the experimental values from 
Ref. [2]. From k 2± 2 x IQ-^SNm^ and G ~ IfcPa, we 
estimate that A ~ 2/im and 1^ ~ 0.3/.tm. The local strain 
in the network induced by a dimensionless displacement 
u is of order u. Thus, for a typical cytoskeletal network 
that becomes nonlinear for strains of order 10-20%, we es- 
timate that r\ is no less than unity, and may be as large 
as 10. Thus, A = ijl^ > 5 since e < 1/2. 



Fig. 5: The transverse displacement u{x) in the buckled state 
of the MT in the presence of higher order nonlinearity from 
numerics. The solid line (black) envelope corresponds to an 
asymptotic exponential (1.4exp(— x/28.8)) fit, with the same 
decay length as for the simple nonlinearity. The material 
parameters are e — 0.01, = !,</) = 3.5. Inset: The circles 
(blue) show u{x) at odd multiples of A/4 (corresponding to the 
peaks in u(x)) as a function of axial length a; on a semilog scale 
for the higher order nonlinearity, while the solid line (black) 
shows the exponential decay with the simple nonlinearity. 



Even for an isolated elastic rod, there are intrinsic 
nonlinearities that arise from the geometry of a buck- 
led shape [17,18]. For small deflections, these nonlinear- 
ities can be obtained by expanding the longitudinal field 
v{x) = ^-^/l -I- (u')2 — 1^ The leading correction to 

what we have considered above is a term of order (u')^, 
which appears in the first and next-to-last terms in Eq.(l). 
Within the limits that we consider here (e < 1, A > 1), the 
dominant effect of these corrections is an additional term 
of order u'^ in Eq.(3). In detail, this geometric nonhnearity 
can be captured within our asymptotic analysis by a sim- 
ple substitution r] —>■ r] + 4>c/'^ near the transition, which 
tends to decrease the amplitude -0 and increase the decay 
length A. But, it is interesting to note that the existence 
of a decay length for the buckling amplitude still depends 
crucially on the existence of finite longitudinal coupling 

"II- 

We find that extrinsic nonlinearities of the em- 
bedding medium can result in a continuous force- 
extension/compression relation, much as intrinsic geomet- 
ric nonlinearities can do for isolated elastic rods [17,18]. 
But, within this model only longitudinal coupling of an 
elastic rod to a surrounding medium can account for a de- 
caying buckling amplitude, such as observed in the exper- 
iments of Ref. [2] . This work also shows how force trans- 
mission by MTs can be controlled by the direct coupling 
of the MTs to the rest of the cytoskeleton, which may be 
regulated by microtubule associating proteins. This can 
allow an embedded MT to bear and transmit supercritical 
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forces over a long range within the cell. Although we have 
considered here only the effects of the confinement of MTs 
by and their coupling to a bulk clastic medium, it is also 
possible that MTs near the plasma membrane may couple 
to the cell substrate. Such coupling can also increase the 
load bearing capacity of MTs. Force transmission such as 
we consider here is different qualitatively from stress prop- 
agation in a homogenous and isotropic clastic continuum 
in that the force is focussed primarily along a single load 
bearing element. We find that this force decays with a 
decay length i to an asymptotic value = /c. For typical 
parameter values of 77 = 1 and e = 0.01, the lengthscale 
£ ~ 35 X ^0 ~ 10/im. Finite compressibility or elastic inho- 
mogeneities in the surrounding medium can also lead to 
finite-range force transmission, although here we consider 
the case of a single stiff filament in an elastic continuum. 
In this case we show that mechanical loads applied at a 
cell's periphery can propagate along MTs throughout the 
cell body without significant diminution. This elucidates 
the importance of MTs in force propagation in a network 
of F-actin with embedded MTs and agrees very well with 
recent experiments [19] that show that presence of MTs 
can significantly enhance the elastic response of such a 
network. While our study focusses on the mechanical re- 
sponse of intracellular microtubules, our results are appli- 
cable to composite elastic media with a wide separation 
of scales in stiffness of its constituents such as biopoly- 
mer networks with bundles and even nanotube-polymer 
composites. 
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